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The one dimensional Kardar-Parisi-Zhang universality class is believed to describe many types 
of evolving interfaces which also have the same characteristic scaling exponents. These exponents 
lead to a natural renormalization group action on the evolution of such interfaces. In this Letter we 
introduce and describe the renormalization group fixed point of the Kardar-Parisi-Zhang universality 
class in terms of a random nonlinear semigroup with stationary independent increments, and via a 
variational formula. Furthermore, we compute the exact transition probabilities using replica Bethe 
ansatz. The semigroup is constructed from the Airy sheet, a four parameter space-time field which is 
the Airy2 process in each of its two spatial coordinates. Minimizing paths through this field describe 
the renormalization group fixed point of directed polymers in a random potential. 
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Interfaces evolving according to local stochastic growth 
rules have been extensively studied in physics, biology, 
material science, and engineering. There has been signif- 
icant theoretical success over the last twenty- five years in 
describing such rough self-affine interfaces which evolve 
due to local processes. In short times such interfaces 
display properties particular to their local processes. 
However, in the long time limit, theoreticians predicted 
that certain universal scaling exponents and exact sta- 
tistical distributions would accurately describe the fixed 
long time properties of a wide variety of rough inter- 
faces. These predictions have been repeatedly confirmed 
through Monte-Carlo simulation as well as experiments. 
What is lacking, however, is a theoretical understanding 
and prediction of the temporal evolution of these inter- 
faces, in the large-time scaling limit. In this Letter we 
describe, both phenomenalogically and statistically, the 
universal temporal evolution of such interfaces. We also 
describe both experimental and numerical means to ver- 
ify these predictions. 

The 1-1-1 dimensional KPZ universality class includes 
a wide variety of forms of stochastic interface growth ^ 
[sl on a one dimensional substrate, randomly stirred one 
dimensional fluids (the stochastic Burgers equation) 
polymer chains directed in one dimension and fluctuatin 
transversally in the other due to a random potential 
and various lattice models such as the driven lattice gas 
model of ASEP and ground state polymer model of last 
passage percolation [y,01- All models can be transformed 
to a kinetically roughening, growing interface reflecting 
the competition between growth in a direction normal 
to the surface, a surface tension smoothing force, and 
a stochastic term which tends to roughen the interface. 
Numerical simulations along with some theoretical re- 
sults have conflrmed that in the long time t scaling limit, 



fluctuations in the height of such evolving interfaces scale 
like t^^^ and display non-trivially spatially correlated in 
the scale t^^^ [J, 0, S| ■ These scales were confirmed ex- 
perimentally in studied involving paper wetting, burning 
fronts, bacterial colonies and liquid crystals [ij. 

Beyond the KPZ scalings, the universality class is char- 
acterized in terms of the long-time limits of the proba- 
bility distribution of fluctuations. These depend on the 
initial data or geometry. Starting from (i) narrow wedge, 
or droplet, one sees the Fque distribution of the Gaussian 
Unitary Ensemble of random matrix theory while start- 
ing from (ii) flat substrate, the -Fgoe distribution of the 
Gaussian Orthogonal Ensemble. A recent series of spec- 
tacular experiments involving turbulent liquid crystals Q 
have been able to not only confirm the predicted scaling 
laws but also the statistics (skewness and kurtosis) for 
the distribution of these fluctuations. The experiment 
involves studying the interface between two topologically 
different turbulent states in a nematic liquid crystal after 
an external voltage (in the form of a laser) is applied so 
as to trigger the Carr-Helfrich instability. 

The multi-point joint distributions of scaled fluctua- 
tions are likewise given by the (i) Airy2 flo| and (ii) Airyi 
11 1 processes. Q showed that the two-point correlation 
function agrees closely with that of the Airy 2 process. A 
further natural initial geometry is two sided Brownian 
motion for which one sees, at later time, a new (though 
correlated) Brownian motion, with a global height shift 
given by the Fg distribution [IJ] . Note that all these spa- 
tial processes have n-point distributions given by Fred- 
holm determinants. 

In this work we answer two questions: 1. Can one 
compute the exact statistics and multi-point joint distri- 
butions for growth off of more general initial geometries; 
and 2. Can one predict multi-time statistics and distri- 
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butions? The answers are provided through a detailed 
study of the KPZ renormalization fixed point which we 
denote as f). We describe the fixed point in two com- 
plementary ways: 1. Through a variational formulation 
similar to that of a stochastically forced Burgers equa- 
tion, but with a new, nontrivial driving noise called the 
Airy sheet; and 2. Through an exact formula for the 
transition probabilities of this Markov process, derived 
by emplo ying the methods of replica Bethe ansatz 13l.ll4| 
(see also [lal). The transition probabilities enable us to 
compute statistics for general initial geometries, and us- 
ing the Markov property, this enables us to also com- 
pute multi-time statistics, thus answering both questions. 
Furthermore, we describe specific experimental and nu- 
merical means to confirm the theoretical results. 

The renormalization operator. Kardar-Parisi- 
Zhang 8] proposed the model equation. 



dth ^ l{d^,h)^ + ^d^,h + ^ 



(1) 



from which the universality class takes its name. 
The noise ^ is Gaussian space-time white noise 
{iit,x)i{s,y)} = 6{t ~ s)S{y - x). Many of the dis- 
crete models have a tunable asymmetry and ^ appears 
as a continuum limit in the diffusive time scale as this 
parameter is critically tuned close to zero [iMl. Re- 
cent progress has been recorded in solving this equation 
[l,,18, 19). © is linearized through the Hopf-Cole trans- 
formation Z — e'^, becoming a stochastic heat equation 
with multiplicative noise dtZ — ^d^Z + ^Z . The narrow 
wedge initial conditions correspond to starting Z with a 
delta function, and the flat initial conditions correspond 
to starting with Z = 1. Correspondingly, in the liquid 
crystal experiments the laser excites a single point, or a 
line. 

The fluctuations are observed under the rescaling, 
h,{t, x) = {Reh){t, x) = ei/2/j(g-3/2^^ ^-i^^^ ^2) 
which satisfies ([1} with renormalized coefficients. 



(3) 



The KPZ fixed point f) is the e — ;> limit of the properly 
centered process /i^. We now construct it. 

The Airy sheet. Let h{u, y; t, x) be the solution of ([1]) 
for times t > u started at time u with a delta function at 
?/, all using the same noise. To center, set h[u,y]t,x) ~ 
h{u, y; t, x) — — log ^2-K(t — u) and define A\ by 

h{u,y;t,x) = + Ai{u,y;t,x). 

After the rescaling ([2]), 

Reh{u,y;t,x) = -^^^ + Aeiu,y;t,x) 

where Ac — ReAi. As e — > 0, Ac{u,y;t,x) converges to 
the Airy sheet A( u, y \ t, x). In each spatial variable it is 
an Airy2 process [lO[. It has several nice properties: 



1. Independent increments. A{u,y;t,x) is indepen- 
dent of A{u', y; t', x) if (u, t) (u' , t') = 0; 
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2. Space and time stationarity. A{u,y; u + 
h,y;t + h, x) '^=^ A{u, y + z;t,x + z); 

3. ScaUng. ^(0, y; i, a;) ^^/^^^(O, t-^/Sy. ^-2/3^.). 

4. Semi-group property. For u < s < t, 

A(ii ir-t T^i — ciim f (x~yf _ {^-yf _ i^-^)'^ 

./il^U, (/, t,XJ — SUPj,gK 2(t-„) 2(5— u) 2(t-s) 

+A{u,y]s,z) + A{s,z]t,x)Y (4) 



KPZ fixed point. Using A{u^ y; t, x) we construct 
the KPZ fixed point t){t,x). By the Hopf-Cole transfor- 
mation and the linearity of the stochastic heat equation, 
the rescaled solution of (H]) with initial data /i° is 



R,h{t 



dy. 



If we choose initial data so that converges to a 

fixed function / in the limit, we can use Laplace's method 
to evaluate t}{t,x) — lim^^Q R^hit^x) ~ To,t/(x) where 

Tu,tfix) := supy {-^^ + Aiu,y;t,x) + fiy)} . (5) 

The operators r„.f, < u < t form a semi-group, i.e. 
Tu.t — Tu,sTs.t which is stationary with independent in- 
crements and T^ t R^X/^Tq iR^-^/^. 

By the Markov property, the joint distribution of the 
marginal spatial process of (] (for initial data /) at a set 
of times ti < <2 < • • • < is given by 

([)(ti), . . . , f)(i„)) {To,tJ, . . . , Tt_„t„ • • • n.tj). 

The process of randomly evolving functions can be 
thought of as a high dimensional analogue of Brownian 
motion (with state space Brownian motions!), and the 
Tti.ti+i as analogous to the independent increments. 

Polymer fixed point. The solution of ([1]) corre- 
sponds to the free energy of a directed random polymer 
x(s), u < s < t starting at y and ending at x, with 
quenched random energy 



H{x{:))^jl{\x{sr-i{s,x{s))}ds 



(6) 



Directed polymers model domain walls in disordered me- 
dia and arise as competition interfaces in multi-species 
growth [20j . Under the rescaling ([2]) this probability mea- 
sure on paths converges to the polymer fixed point; a con- 
tinuous path ■nu,y;t,x{s)^ u < s < t from y to x which at 
discrete times u = sq < ■ ■ ■ Sm-i < t is given by the 
argmax over xq, . . . , a^m-i of 

{Tu,siSy){xi) + {Ts^,s2Sxt){x2) H h (Ts,„_i,t4„,_i)(a:)- 

(7) 

This is the analogue in the present context of the mini- 
mization of the action and the polymer fixed point paths 
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FIG. 1: Polynuclear growth with two nucleations distance e ^ 
apart observed at times e~^^^ti {i = 1,2,3) in an spatial 
scale and e"^''^ fluctuation scale. 



are analogous to characteristics in the randomly forced 
Burger's equation. One might hope to take the analogy 
farther and find a limit of the renormalizations of ([5]), 
and minimize it to find that path T^u.y-t.x- However, the 
limit does not appear to exist, so one has to be satisfied 
with the limiting paths themselves. The path no^y-t^x 
turns out to be Holder continuous with exponent 1/3—, 
as compared to Brownian motion where the Holder ex- 
ponent is 1/2—. As the mesh of times is made finer, a 
limit £{TTo^y-t,x) of (O does exist, and through it we can 
write the time evolution of the KPZ fixed point in terms 
of the polymer fixed point through the analogue of the 
Lax-Oleinik variational formula, 

i){t,x) = SUPyf,jf{£{7ro,y;t,x) + f{'^0,y;xA^))}- (8) 

Universality. The KPZ fixed point, Airy sheet, and 
polymer fixed point are universal and will arise in random 
polymers, last passage percolation and growth models 
- anything in the KPZ universality class. Just as for 
([1]), at the microscopic scale, approximate versions of the 
variational problem ([5]) hold, becoming exact as e 0. 

For example, consider the PNG model [13] with a finite 
collection of nucleations spaced order apart (Figure 
[T]). At time e~^/'^t we look at e^^/^ scaled fluctuations 
in spatial locations e~^x. As e goes to zero, these fluc- 
tuations converge to f) where the initial data / is — cxd 
except at the nucleation points, where it is zero. By in- 
troducing additional nucleations at time on the order of 
e^^/^ and spatial locations order apart, it is possible 
to modulate the value of / at these non — oo points. Tak- 
ing the number of nucleation points large allows one to 
recover any /. This scheme provides a easy to implement 
approach to run numerical simulations of the fixed point 
for general initial data. 

The experiment of Q is well described by the KPZ 
fixed point with a single nucleation. Future experiments 
could probe the effect of additional nucleations. Using 
statistics to differentiate between types of initial data 
given finite time observations is a driving force for the 
development of the following exact formulas which pro- 
vide theoretical predictions. 



Transition probabilities. We start with a simple 
version. Given Ci,yi G M, i = 1, . . . ,i>, and Sj, Xj e M, 
J = 1,...,M, let l){t,x) = {Ttf){x) with f{y,) = -c, + 
iVi — xi)'^ /2t and / = —oo otherwise. Then, 



P(f)(i,a;,)<s,- 



,/Lt) = det(/ - L) 



where L is the operator with kernel 

L(z,z')= / dvi---dv^+idui---dUf,e^y^~''''>"K\vi) 

Ja{s,c) 

acting on L^(]R), where H = —j^ + x is the Airy oper- 
ator, K is the Airy2 operator which is the spectral pro- 
jection of H onto its negative eigenvalues, and 

v4(s, c) = {(u, v) : max{ui - s;} -(-max{wj — c^} > Vu+i}- 

1=1 j=i 

We assume xi < X2 < ■ ■ ■ < Xf^ and yi > ■ ■ ■ > Vu- The 
formulas give a consistent family of finite dimensional 
distributions. Specializing to ^ — \ oy v — \ one checks 
that the resulting process (in xi and yi) is the Airy2 
process in each variable. 

From the above the continuous transition probabilities 
can be obtained by approximation. Let Y^^ ^ be the oper- 
ator defined by Yj^ b\f{y) — '^{^i v) where u is the solutio 
of the boundary value problem, dtu — —Hu for a < t < b, 
with Dirichlet condition u(t, x) = on x > g(t) and ini- 
tial data u{a,x) = f{x). 

Let c{x) and s{x) be finite in [£c,rc\ and [£s,?'s] and 
-l-oo outside. Let f{x) — — [c{x) — {x — £s)'^/2t^, g{x) = 
s{x) - {x~lsf/2t and 

p{0J;t,9) = P{To,tf{x) < g{x), - oo < a; < oo) 

Then p{0, f;t, g) is given by a Fredholm determinant 
det(/ — K + L) where L{z, z') is given by 



{z[u) 



where c{y) — —f{£c + fc — y), and K, H are as before. 

Although the connection is not yet understood, our 
formula should be a scaling limit of a less explicit de- 
terminental formula derived earlier for the microscopic 
model TASEP [ill, 

Replica Bethe ansatz. Given Sj,Xj, j — 
and Ci,2/i, i = 1, . . . , consider J2j=i Z{t^Xj) where 
Z{t,x) solves the stochastic heat equation with initial 
data X]r=i ^~'^'^yi- Fi'om the linearity this can be writ- 
ten as 'Ek,i=i<^~'"'~'" Z{yhXk;t). Following [H [l^ we 



compute the generating function 



G{s,x;c,y) 



exp 



(-e24 J2 



'''--^Z{yi,Xk;t)\)) 



By studying asymptotics as t goes to infinity we obtain 
the formula for the transition probabilities. 
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Expanding the generating function exponential we 
write this using rephcas as 



G{s,x]c,y) = 1 + Ew=i — -W. 



ATI 



N\ 



where the N particle wave functions of the (5-Bose 
gas Hn = -\Y.f=idl, - \ Y.^^j=i^{x^ - Xj) are 

U^l =Er„..i„=ie"'^'^+-+^'"^(y^i---2/^J and = 

■■■xi^). The wave func- 
tions are symmetric, so the propagator is only needed 
on the symmetric subspace. Thus we may employ the 
eigenfunction expansion of iJ„ ^13i] . Expanding the ten- 
sor products given above, 



G{s,x]c,y) = 1 



N=l 



(-1)' 



AT! 



y-M „-(^fci+---+^fc,„) ^'-(^fcl---'^'=Ar) 
Z^fel,...,fcjv = l 

X^" p-(cii+---+ci„) 

Z^;i,...,/jv=i 



jlriVli ,---,Vlj. 



We use the factorization assumption of jl4| . 
Z^fei,...,/cjv=i ^ Vv^o) 



(9) 
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which holds asymptotically in t (as confirmed by the re- 
covery of the Airy2 process in the limit when v = 1 [l^ ) . 
From this we obtain an asymptotically correct generat- 
ing function given by a Fredholm determinant det(l — N) 
where N{z^ z') is given by 



/2t^y-^ /2t'^z,z'>oe-'y e-^^jduM{u+z)k\{u+z')^{ 



au — s:-c 



where for vectors s = (si, . . . , s^) and c = (ci, . . . , c^), 

Cj - i+(esi+...+eSj,)(eci+...+eo^) ^^^^ - 

^ X]a=i ^a^sa • Jy is the analogue with s, x replaced 
by c, -y, and T^2/2t,y2/2t = exp{^ X^Li ^^a^^a + 
it^ ^1=1 Va^ca} is the parabolic shift. After a series 
of integrations by parts JJ| we arrive at 

Conclusions. We have described new random ob- 
jects: A Markov process, the KPZ fixed point described 
through a variational formula on a network of minimiz- 
ing paths, the polymer fixed point. Both are described 
in terms of the Airy sheet. Our description implies vari- 
ational relations between the Airy processes, which are 
known to appear as spatial interface fluctuation processes 
for this class, and extend the universal limit processes of 
random matrices. 

Furthermore, we have computed exact transition prob- 
abilities for the KPZ fixed point in terms of Fredholm 
determinants and prescribed both experimental and nu- 
merical methods to verify our prediction. For exam- 
ple, exciting several nucleations in turbulent liquid crys- 
tals and studying the correlations of interface height 



fluctuations at different points. These fixed points are 
expected to describe the asymptotic fluctuations over 
large space and time of a very wide class of one di- 
mensional physical systems with nonlinear transport of 
fluctuations, such as those modeled by directed random 
polymers, random growth models, stochastic Hamilton- 
Jacobi and Hamilton- Jacobi-Bellman equations, as well 
as the Kardar-Parisi-Zhang equation. 
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